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ABSTRACT 

A search for young massive star clusters (YMCs) in the nearby face-on spiral galaxy 
M51 (NGC 5194) has been carried out using UBV CCD images from the prime focus 
camera on the Lick 3 meter Shane telescope. The YMC population is found to be 
quite rich with a specific J7-band luminosity Tl[U) ^ 1.4, consistent with the high 
current star formation rate of this galaxy. The brightest clusters have My ~ —12.5, 
far brighter than any young clusters currently known in the Milky Way and even 
surpassing the luminosity of the R136 cluster in the 30 Dor complex in the Large 
Magellanic Cloud. A few of the YMCs are examined on archive HST/WFPC2 images, 
confirming their cluster nature and providing estimates of their effective radii of 2 — 3 
pc. The number of YMCs in M51 is compatible with extrapolation of a power-law 
luminosity function with exponent ~ — 2 from a Milky Way- like population of open 
clusters. Both the SFR and Tl(?7) of M51 are similar to those of other cluster-rich 
spiral galaxies hke NGC 1313 and M83. 
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1 INTRODUCTION 

M51 is probably among the most famous galaxies in the sky, 
instantly recognizable by its nearby companion NGC 5195. 
Also known as the "Whirlpool Galaxy", it was one of the 
first galaxies in which spiral structure was discovered by 
Lord Rosse, and it provides a spectacular text-book example 
of a grand-design Sbc type spiral seen nearly face-on. It ha s 
a very high optical surface brightness (Okamura et al. 197C) 
and is also a strong emitter of far-infrared radiation, i ndicat- 
ing s t rong star formation ac tivity in the spiral arms (Smith 
1982; Hippelein et al. 1996). Deep images reveal fain t out- 



lying mate rial embedding both M51 and NGC 5195 (Burk 



head 1978), providing clear evidence that the two galax- 
ies are physically interacting. These characteristics, along 
with the relatively small distance (8.4 ± 0.6 Mpc, Feldmeier 
et al. 1997) and the location far from the galactic plane 
(b = 68°) make M51 an attractive target for studying its 
population of Young Massive Star Clusters (YMCs). 

YMCs are abundant in interacting and merger galax- 
ies like e.g. the "Antennae" NGC 4038/4039 (a list of 
galaxies with known YMC populations is given in Larsen 
(1999b)). However, th ev are also seen in rn ore normal galax- 



ies like e.g. the LMC (jvan den Bergh 199l|) and M33 ( Phris- 



tian & Schommer 198^). In a recent study of 21 n earby . 



non- interacting spiral galaxies, Larsen & Richtler (200C) 
(LR2000) found a strong correlation between the specific 



[/-band luminosity Tl{U) of YMCs in a galaxy and the area- 
normalized star formation rate Esfr- In this respect, M51 
provides a highly interesting intermediate case of a clearly 
interacting galaxy that has still retained the characteristics 
of a normal spiral. 

This paper reports the results of a study of YMCs in 
M51. In Sect. ^, YMCs are identified and photometry is 
obtained from ground-based UBV CCD imaging. Archive 
HST/WFPC2 images are then used to examine a few clus- 
ters in detail, including measurements of their sizes (Sect. 
Next, the specific [/-band luminosity of the M51 cluster sys- 
tem is derived and compared with other galaxies (Sect. ^). In 
Sect, ^the YMC population in M51 is compared with young 
clusters in the Milky Way and the LMC, along with some 
considerations on formation of YMCs. Finally, conclusions 
are in Sect. ^. 



2 DATA 

2.1 Observations and initial reductions 

CCD images in the UBV passbands were obtained on Mar 
13 - 14, 2000 with the prime focus camera (PFCAM) on the 
Lick 3 meter Shane telescope at Mount Hamilton, California. 
The total integration times were 3600, 1200 and 900 sec. 
in U, B and V , split into 3 individual exposures in each 
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Table 1. HST archive images of M51 used for photometric cah- 
bration of ground-based photometry. 

PID PI Filter Exposure times 

5777 Kirschner F555W 600 s 

5777 Kirschner F439W 2 X 700 s 

5652 Kirschner F336W 3 X 400 s 



filter. Typical seeing values ranged between I'.'O and 1'.'5, 
but the image quality was degraded by problems with the 
alignment of the mirror cell after recent re-aluminization 
and oscillations due to wind, making stellar images appear 
somewhat elongated. The image scale was 0'.'296 /pixel and 
the total field of view was 10 X 10 , sufficient to fully cover 
both M51 itself and the companion NGC 5195. 

To reduce read-out time, the two halves of the CCD 
were read out in parallel through two amplifiers. During the 
subsequent reductions, the difference between the bias levels 
in the two parts of the CCD image was found to change by up 
to ^ 40 ADU from one exposure to another. The difference 
was eliminated by adding a constant number to all pixel 
values in one half of the image. A more serious concern was 
that the fiat-field varied significantly, most likely because of 
scattered light. Dividing two skyflats in the same filter with 
each other, the ratio was found to vary by up to ~ 20% from 
the centre to the corner of the chip, with the most severe 
problems in the southern ~ 300 pixels (about l!5) of the 
field. Within the central 1000 x 1000 pixels the variations 
were less dramatic, but still quite significant at about the 
5% level. The effect was clearly visible in the final calibrated 
images as large-scale gradients in the background. However, 
the pattern appeared to be more or less the same in all 
passbands so the fiatfielding errors may, to a certain extent, 
cancel out for colour indices. 



2.2 Photometric calibration 

Standard fields from Landolt ( 1992| ) were observed for pho- 
tometric calibration. The transformation equations were as- 
sumed to be of the form 



NGC 5194 photometry 



V = V + X {b - v) + Zy 
B~V = C(f,_„) X (6 - t)) + 
U~B = C(„_i,) X (u - 6) + 

where capital letters denote standard magnitudes and small 
letters indicate instrumental magnitudes, corrected for at- 
mospheric extinction. However, because of non-photometric 
conditions and the flat-field problems it was necessary to 
correct the zero-points of the UBV photometry using HST 
archive images. UBV magnitudes were obtained for sources 
in the HST images listed in Table 0, followin g the stan- 
dard procedure described in Holtzman et al. ( 1995| ). Ob- 
jects bright enough to be used for calibration of the ground- 
based images were easily measured through the Holtzman et 
al. reference aperture of r = 5 pixels, avoiding problematic 
aperture corrections. The zero-points {zv, 2:(i,_„) and 
of the ground-based magnitudes were then adjusted to fit 
the HST system, while the scaling constants (c„, C(i,_„) and 
C(ti-6)) derived from the Landolt standard fields were pre- 
served. The overlap between the HST datasets used for the 
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Figure 1. {B — V, V) colour-magnitude diagram for objects in 
the M51 field. Objects with B — V> 0.5 are mostly foreground 
stars whereas contamination blueward of B — V = 0.5 is expected 
to be small. 



BV and U calibrations amounted to less than the area of 
one WF camera chip, making it difficult to find a suitable 
number of calibration objects. The U — B calibration is thus 
somewhat more uncertain than the B — V and V^-band cali- 
brations and it is estimated that is good to ±0.15 mag 
while 2:{i,_„) is probably accurate to better than ±0.05 mag. 
Although significantly worse than what could have been ob- 
tained under optimal photometric conditions, the calibra- 
tions adopted here are nevertheless accurate enough for the 
conclusions drawn in this paper. 



3 SELECTION OF CLUSTER CANDIDATES 

In order to facilitate a comparison with the results of 
LR2000, the same data reduction and analysis procedures 
were followed as closely as possible. Obj ects were dete cted 
using the daofind task in daophot ( Stetson 1987| ) on 
background-subtracted V and B band images and match- 
ing the two object lists. Aperture photo metry was then car- 
ried out as described in Larsen (1999a), using an aperture 
radius of r = 5 pixels. Finally, a correction for interstellar 



absor ption oi Ab = 0.150 mag was applied ( gchlogcl et al 
199^ ). 

The resulting {B — V,V) colour-magnitude diagram is 
shown in Fig. 0. An initial list of cluster candidates was 
then obtained by selecting all objects with 



B 

Mv < 



V < 0.45 

-9.5 for U-B < -0.4 
-8.5 for U-B > -0.4 



These criteria ensure minimal contamination from fore- 
ground stars, the majority of which will be redder than the 
B — V cut, and by individual luminous stars in M51 itself 
because of the Mv limit. However, many of the objects iden- 
tified in this way were still not likely cluster candidates, but 
high-surface brightness parts of spiral arms, HII regions or 
just artefacts found by daofind where no real object ex- 
isted. A visual inspection of the images, looking for rela- 
tively isolated, point-like sources, was found to be the only 
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Figure 2. The intrinsic FWHM vs. U — B colour for cluster 
candidates in M51 that passed the initial visual inspection. The 
horizontal dashed line indicates the size cut applied to further 
constrain the cluster sample. 



way to isolate good cluster candidates from the multitude 
of other objects in the frames. 

LR2000 further used Ha images to reject the youngest 
objects, still embedded in HII regions and often surrounded 
by a fuzz of nebular emission and star forming regions, but 
no Ha filter was available at PFCAM during this observ- 
ing run. Instead, such objects were rejected based on their 
angular siz es. Object siz es were measured using the ishape 
algorithm (Larsen 1999a), convolving the point spread func- 



tion (PSF) of the 1/— band image with King profiles to match 
the observed object profiles. A concentration parameter of 
30 was assumed for the King profiles. Fig. ^ shows the in- 
trinsic object sizes derived by ishape as a function of U — B 
colour for objects that passed the initial visual inspection. 
Clearly, most of the more extended sources are found among 
the bluest objects, while most objects with U—B ^ —0.6 
have FWHM < 2 pixels. Note that the typical seeing in 
our images corresponds to FWHM ~ 5 pixels. A size cut at 
FWHM=3 pixels was finally applied to exclude the fuzziest 
objects, likely to be OB associations or HH regions rather 
than real star clusters. 

Although the cluster colours have been corrected for 
Galactic foreground extinction, some residual reddening 
may still be present from M51 itself. The question of the 
optical thickness of face-on spiral galaxies is still debated, 
and for any particular object the absorption will obviously 
depend on the local structure of the s urrounding interstel - 
lar medium. However, a recent study (Xilouris et al. 199E) 



concluded that face-on spiral galaxies are essentially trans- 
parent with an optical depth of less than one in all opti- 
cal bands. M51 being obviously face-on (inclination angle 
~ 20°, TuUy 1974), interstellar extinction problems are thus 
expected to be small. 

The final list consists of 69 cluster candidates (Table |^. 
A {B — V, U — B) diagram for these is shown in Fig. ^, along 
with stellar mode ls by Bertelli et al. (1994) and the Gi- 



rardi et al. ( 1995 ) "S" -sequence which marks the average 
colours of LMC clusters. The arrow indicates the effect of a 
reddening of 0.3 in B — V, equivalent to 0.93 magnitudes of 
y-band absorption for a standard reddening law. The M51 
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Figure 3. {B — V , U — B) two-colour diagram for the final sample 
of cluster candidates in M51 



clusters tend to have somewhat bluer U — B colours than ex- 
pected from the S-sequence, but a shift of about 0.1 in U — B 
would make them match the S-sequence colours nicely. Such 
a shift would be consistent with the estimated uncertainties. 
The S -seque nce is essentially an age sequence and Girardi 
et al. (1995) published a calibration of age as a function of 
-B) diagram. However, because of 



position in the {B — V ,U - 
the uncertainties in the photometric calibration and espe- 
cially in the U band, a detailed discussion of ages for the 
M51 clusters will not be given in this paper, but note that 
the clusters do scatter along the S-sequence, indicating a 
significant age spread. Because of the fading with age, the 
sample becomes increasingly incomplete for higher ages. The 
upper age limit, due to the B — V cut, is about 500 Myr. 

Fig. ^ shows a V-band image with the cluster positions 
indicated. Most of the cluster candidates are located in or 
near the spiral arms of M51 itself. This is to a large extent 
a selection effect, due to the fact that the youngest clusters 
are more luminous (for any given mass) and therefore easier 
to detect. However, note that also the companion galaxy 
contains a couple of massive clusters. The HST F439W + 
F555W pointing is also shown in Fig. ^ 



4 HST DATA 

In addition to being useful for calibrating the photometry, 
the superior angular resolution of the HST data provides 
a welcome way of examining the structure of the cluster 
candidates in much greater detail than what is possible from 
the ground. Fig. ^ shows WFPC/2 close-ups of 10 clusters, 
marked by asterisks (*) in Table ^ The images in Fig. ^ are 
from two combined F439W exposures, allowing for cosmic 
ray rejection. 

First, the HST images clearly confirm the cluster nature 
of most of the objects. Exceptions may be objects #627 - 
a quite faint source located near the nucleus of M51, and 
#597 and #403 which may consist of two or more closely 
separated individual stars. The remaining seven objects are 
far too bright to be individual stars and appear well resolved. 

Most of the clusters are evidently not isolated objects, 
but are surrounded by individual luminous stars and other 
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Table 2. Comparison of ground-based and HST- based data for 7 clusters in M51. Cols. 2 and 3 list ground-based V and B — V magnitudes, 
Cols. 4 and 5 give the corresponding values measured on HST WFPC/2 images and the differences AV = V{ground) ~ V(HST) and 
AB — V = S — y(ground) — B — V(11ST) are in Cols. 6 and 7. The colours given here have not been corrected for reddening. Effective 
cluster radii measured on the HST images are given in the last column. 



ID 


V (ground^ 


B — V (^ground) 


V (HST) 


B — V fHRT"! 


Ay 


AB — V 




403 


19.58 


0.342 


19.80 


0.29 


-0.22 


0.05 


2.7 


477 


18.41 


0.210 


18.46 


0.33 


-0.05 


-0.12 


0.9 


479 


18.11 


0.098 


18.83 


0.06 


-0.72 


0.04 


3.8 


534 


19.48 


0.177 


20.24 


-0.04 


-0.76 


0.22 


2.3 


589 


19.40 


0.070 


19.51 


0.08 


-0.11 


-0.01 


2.1 


597 


19.90 


0.177 


19.72 


0.52 


0.18 


-0.34 




627 


20.13 


0.237 


20.53 


0.51 


-0.40 


-0.27 


0.0 


693 


19.42 


0.148 


19.66 


0.11 


-0.24 


0.04 


2.7 


7f3 


18.33 


0.071 


19.22 


-0.11 


-0.89 


0.18 


0.8 


9f2 


19.78 


0.317 


20.34 


0.16 


-0.56 


0.16 


2.2 




403 



477 



479 



534 



Figure 4. A V— band image of M51 with the cluster positions in- 
dicated. South is up. The large-scale variations in the background 
are due to poor flatfielding, presumably because of scattered light 
in the skyflats. The F439W / F555W HST pointing (PID 5777) 
is also shown. 



clusters. This is, however, to be expected because most of 
these clusters have very blue U — B colours, indicating low 
ages. A U — B colour of —0.8, typical for the clusters in 
Fig. 1^, corresponds to an age of about 12 Myr (Sirardi et 
al. 199^), too short for a young star cluster to escape from 
the environment where it was born. The two reddest clusters 
(#693 and #912) with U-B = -0.59 and U-B = -0.58 
aheady appear somewhat more isolated than the other clus- 
ters in the Fig. ^ consistent with the higher ages inferred 
from their U — B colours (~ 30 Myr). 



4.1 HST versus ground-based photometry 

Table ^ lists HST-based photometry and sizes (see Sect. 4.2) 
for the 10 objects in Fig. ^ For comparison, ground-based 
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Figure 5. Clusters identified on HST images. The circles are 2 
arseconds in diameter. 



photometry is also listed. The photometry in Table |2| has 
not been reddening corrected. 

Magnitudes measured on the ground-based images are 
generally brighter than those based on HST images. This is 
not surprising, considering that the HST magnitudes were 
measured through a 0'.'5 aperture which clearly does not in- 
clude all the hght surrounding the clusters that will enter 
into a ground-based aperture. Generally, the agreement on 
B — V colours is better than on V magnitudes, confirming 
that colours can be more accurately measured than individ- 
ual magnitudes for objects in crowded fields. The difference 
between ground-based and HST-based B — V colours clearly 
increases for the fainter objects, with the largest discrepancy 
for the likely "non-clusters" (#597 and #627) pointed out 
above. 

It is also worth noting that the HST pointings were 
centred on the nucleus of M51, providing a "worst-case" 
comparison between HST and ground-based data. Further 
out in the disk of M51 the surface brightness drops and the 
cluster sample becomes less biased towards the brightest 
(and thus youngest) objects where crowding problems are 
most significant. 
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Figure 6. Histogram of the size distribution for objects in tiic 
WF4 chip of the combined F439W image. 



4.2 Cluster sizes 

Sizes were measured on the WFPC/2 F439 W imaRes for th e 
clusters in Fig. ^ with the ishape algorithm (Larsen 1999a). 
At the distance of M51, one WF pixel corresponds to a linear 
scale of 4 pc, so measuring cluster sizes of the order of a few 
pc is not trivial. The input PSF for ishape was generated 
empirically from point sources in the HST images, using the 
PSF task in daophot. This was preferred o ver generating 



the P SFs with the TinyTim PSF simulater ( Krist fc Hook 



1997 ) which does not include the so-called "diffusion kernel" 



for subsampled PSFs, and would thus have caused ishape 
to derive too large cluster sizes. 

The clusters were modelled as King profiles with con- 
centration parameter c = 30, resulting in the half-light radii 
{Re) listed in the last column of Table ^ The cluster sizes 
measured here are consistent with sizes for YMCs in merger 
galaxies like the Antennae (Whitmore et al. 1999), with 
_Re ~ 2 — 3 pc and do also match the typical size of globu- 
l ar clusters i n the Milky Way and other galaxies quite well 
(Harris 199C). It has been estimated from simulations that 



ISHAPE is able to measure sizes down to times the size 



of the PSF when sufficient signal is present (Larsen 1999a), 
corresponding to ~ 0.5 pc in the current case. To illustrate 
the reliability of the cluster sizes. Fig. ^ shows a histogram 
of the size distribution for all objects brighter than V = 22 
in the WF4 chip. The histogram shows a narrow peak at 
-Ro = 0, corresponding to unresolved point sources, while 
objects with sizes characteristic of the YMCs appear well 
separated from the -Re = peak. The YMCs are thus well 
resolved and their sizes are well determined. 



5 GLOBAL PROPERTIES OF THE M51 
CLUSTER SYSTEM 

5.1 Specific luminosity 

In studies of old globular cluster systems (GCSs) it has be- 
come customary to characterize the richness of a GCS by 
its spec ific fr equency, originally defined by Harris & van den 
Bergh (1981) as the number of globular clusters (Nqc) per 
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Figure 7. U — B vs. B — V for galaxies with T-type in the range 
3 to 5 (from RC3). The dashed line is a fit to the data points. 



host galaxy luminosity (Mv): Sn = Nqc x io0.4x{Mv+15) ^ 
Howev er, a potenti ally better measure is the specific lumi- 
nosity (Harris 1991), which is the ratio of the total luminos- 



ity of the cluster system to that of the host galaxy. Specific 
luminosities are much less sensitive to incompleteness effects 
at the lower end of the cluster luminosity function because 
the fainter clusters obviously contribute less by light than 
they do by number. For YMC systems this advantage be- 
comes even greater than for old GCSs because the faint wing 
of the cluster luminosity function is usually unknown, so a 
definition of an equivalent specific frequency for YMCs has 
to rely on artificial and somewhat arbitrary magnitude lim- 
its. 

LR2000 defined the specific f/-band luminosity for 
YMCs as follows: 



Tl{U) = 100 . 



i'Galaxy(C^) 



(1) 



where I/ciustors((7) and LoaiaxyiU) are the total [/-band lu- 
minosities of the cluster system and of the host galaxy, re- 
spectively. The U band was chosen to sample the young clus- 
ter and stellar populations in a galaxy as cleanly as possible. 

To calculate Tl{U) for M51, Lciustcrs(f/) was obtained 
simply by adding up the [/-band light from all clusters in 
Table |^. Unfortunately, no integrated f/-band photometry 
is available for M51 itself so its [/-band magnitude was esti- 
mated from the BV photom etry given in the RC3 catalogue 
( |de Vaucoulcurs et al. 1991 ). Fig. |^ is a plot of reddening- 
corrected {U — B)o vs. {B~V)o for all galaxies in RC3 with 
T types between 3 and 5, M51 itself being classified as Sbc 
(T = 4). The data in Fig. can be fitted by a straight line, 
yielding 



{U-B)o = 1.01 X {B-V)o - 0.59 



(2) 



with a scatter of a{U — B) = 0.09 around the fit. The in- 
tegrated {B — V)o colour for M51 is given as 0.53 in RC3, 
so from Eq. (|) we get {U-B)o = -0.05. This should 
be accurate to 0.1 magnitudes, translating into a 10% un- 
certainty on the derived Tl{U) value. These {B — V)o and 
{U — B)o colours also compare well with integrated photome- 
try for NGC 5236 (M83), a galaxy which is in many respects 
quite similar to M51. For NGC 5236, RC3 gives integrated 
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{B — V)o and {U — B)o colours of 0.61 and —0.01, respec- 
tively, i.e. M51 is slightly bluer in both B~V and U~B. 

The specific [/-band luminosity of the M51 YMCs can 
now be obtained: 



Tl{U) = 1.36 ±0.28 



(3) 



Similarly, the V^-band specific luminosity is found to be 
Tl{V) = 0.41 ± 0.07. The errors quoted here are the for- 
mal errors due to sample statistics. The real uncertainties 
are probably higher, but are hard to estimate. For exam- 
ple, lowering the size cut to FWHM=2 pixels instead of 3 
pixels would change Tl{U) to 0.95 ± 0.24. On the other 
hand, the M51 photometry does not go quite as deep as 
most of of the galaxies in the sample of LR2000, - primar- 
ily because M51 is further away. Thus TLiU) and Tl(F) 
are probably both somewhat underestimated. Tl{V) may 
be more affected because older, redder clusters carry more 
weight than for Tl{U), while at the same time being fainter 
and therefore more incomplete. 
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Figure 8. Specific [/-band luminosity as a function of star for- 
mation rate surface density for galaxies in the sample of Larsen 
& Richtler (2000) and M51 (NGC 5194). 



5.2 The TLiU) vs. Esfr correlation 

LR2000 found a very good correlation between Tl{U) and 
the area-normalized star formation rate Esfr for the 21 
galaxies in their sample and an additional 10 galaxies 
(mostly starbursts and merger galaxies) for which data were 
collected from the literature. How does M51 fit into this re- 
lation? 

To derive Esfr for M51, we use the relation from 
LR2000: 

EsFR(Moyr-'kpc"') = 144000 x iQ-O-^mFm - 2 log Do ^ 

based on the calibration in Buat & Xu ( |l996| ). rriFiR is the 
RC3 FIR magnitude (based on IRAS 60^i and 100^ flux 
densities) and log Do is the logarithm of the optical galaxy 
diameter in tenths of arcminutes. Note that both Esfr and 
Tl{U) are distance independent. In some cases the RC3 "ef- 
fective aperture" (logAe), containing half the _B-band light, 
may provide a more reasonable measure of the area to which 
the SFR should be normalized, but this is not catalogued for 
all galaxies in RC3 and thus log Do was chosen for homo- 
geneity. 

For M51, mFiR = 7.86 and log Do = 2.05. Inserting 
these numbers in (^, the star formation rate surface density 
is then 



Esfr = 8.21 x lO'^M© yr"^ kpc" 



(5) 



This number is rou ghly a factor of two lower than the value 
given in Kennicutt (19981), derived from the Ha luminosity. 



Considering the intrinsic uncertainties in both methods, the 
agreement is quite satisfactory. 

The Tl{U) - Esfr plot in LR2000 can now be updated 
with one more data point for M51, as shown in Fig. |^. M51 
(labelled by its NGC number, 5194) fits quite nicely into the 
existing relation, with a Tl{U) comparable to other high 
Esfr galaxies like NGC 1156 and NGC 1313. It may faU 
somewhat below the Tl{U) value expected for its Esfr, but 
the error bars in the horizontal direction of Fig. ^ are prob- 
ably substantial, with an uncertainty of at least 50% on the 
calibration of star formation rate as a function of FIR lumi- 



6 DISCUSSION 

6.1 Comparison of the luminosity functions of 

YMCs in M51 and open clusters in the Milky 
Way 

M51 contains a large number of very luminous star clusters, 
including 40 clusters with Mv < — 10 and 4 clusters with 
Mv < -12 (Table |). As indicated by Table |, the ground- 
based V magnitudes are systematically brighter than the 
HST photometry, with an average difference of AV — —0.38 
for the objects listed in Table H. Applying a correction of 
0.38 mag to all V magnitudes, M51 still contains 30 clusters 
with Mv < —10 and 3 with V < —12. For comparison, 
the brightest known young clusters in the Milky Way have 
My ~ -10 {h and x Persei, Schmidt-Kaler 1967; NGC 3603, 
van den Bergh 1978). Although it cannot be excluded, of 
course, that a few highly luminous clusters might be hiding 
in remote parts of the Galactic disk, it appears unlikely that 
the Milky Way contains any significant number of highly 
luminous clusters similar to those in M51. The brightest 
clusters in M51 even outshine the R136 cluster at the centre 
of the 30 Dor region in the LMC which has My = -10.6 
(Ivan den Bergh 19781). 



Do YMCs represent a separate class of objects, distinct 
from low-mass open clusters, or can YMC populations be 
reconciled with normal populations of low-mass clusters, as- 
suming a standard power-law luminosity function (LF)? Let 
us compare the YMC population in M51 with the system 
of open clusters in the Milky Way, for the moment ignoring 
other differe nces b etween the two galaxies, van den Bergh & 
Lafontaine ( 1984 ) studied the luminosity function of Milky 
Way open clusters and found it to be well represented by 
the relation 



logcr(Afy) = 1.55-^0.2My 



(6) 



nosity alone ( Buat fc Xu 199!: ) 



for clusters fainter than Mv = —9, where a(Mv) is the 
number of clusters per kpc^ per magnitude bin. In luminos- 
ity units, this becomes 

n(L) (X (7) 
where L is the luminosity, n(L) is the number of clusters per 
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luminosity bin and a = —1.5 (Elmegreen & Efremovl997). 
In the Milky Way, extrapolation of to higher luminosi- 
ties would yield a total of ~ 100 clusters with Mv ~ —11, 
assuming that the Galactic disk has an area of 500 kpc^. 
This is clearly incompatible with wh at is a ctually observed 
and van den Bergh & Lafontaine (1984) suggested that 
the luminosity function of open clusters in the Milky Way 



Even higher levels of star form ing activity are seen in 
merger galaxies l ike the Antennae (Whitmore et al. 199£) 



— 11 < Mv < —8. In fact, there must be a drop-off in a lu- 
minosity function of the form (^) at some magnitude level, 
as the total luminosity of a cluster system with a power-law 
LP diverges for a > —2. However, it is less obvious that the 
drop-off will occur at the same magnitude in all galaxies. 

Young LMC clusters seem to follow a similar power- 
law LF with a — —1.5 ± 0.2, but the LF of young clusters 
declin es more slowly at br ight magnitudes than in the Milky 
Way dElson fc Fall 198£| ) due to the LMC's population of 



YMCs. So at least in the LMC, YMCs appear to mark a 
natural extension of the open cluster population to brighter 
magnitudes. 

M51 contains 23 and 13 clusters in the two magnitude 
intervals -11 < My < -10 and -12 < My < -11, respec- 
tively. If the AV = 0.38 correction is applied, the numbers 
are 20 clusters with — 11 < My < — 10 and 6 clusters with 
-12 < Mv < -11. Most of the visible part of the disk of 
M51 is confined within a diameter of about 6' correspond- 
ing to an area of 170 kpc^. Applying Eq. (h) to the two 
magnitude intervals above, 60 and 38 clusters would be ex- 
pected - three times more than what is actually observed. 
Of course, the extrapolation of Eq. to such bright magni- 
tudes is very sensitive to the exact value of the exponent a. 
Changing a to —2 instead of —1.5 and using the Milky Way 
clusters in the magnitude bin — 5 < Afv < — 8 for reference 
(5.2 clusters kpc~^, van den Bergh and Lafontaine 1984), the 
predicted numbers drop to 6 clusters with — 11 < My < — 10 
and 2 in the interval —12 < My < —11, now several times 
below the observed numbers. 

Although very crude, these calculations show that even 
though M51 is apparently much richer in luminous clusters 
than the Milky Way, the number of low-luminosity clusters 
in the two galaxies could still be quite similar without requir- 
ing any particularly peculiar luminosity function. The YMCs 
in M51 could easily be accomodated by simple extrapolation 
of a power-law luminosity function with an exponent some- 
where between —1.5 and —2, and there is apparently no need 
to conceive them as a separate class of objects. 



6.2 Formation of massive clusters 

From Fig. |^, M51 is indeed quite rich in YMCs compared to 
the "average" spiral galaxy. However, other spiral galaxies 
with similar rich YMC populations are known, and there is 
no particular evidence that the presence of YMCs in M51 
is boosted by the interaction with NGC 5195. If anything, 
M51 may be slightly cluster-poor for its Esfr,. None of 
the other cluster-rich galaxies in the LR2000 sample show 
obvious signs of interaction, and the differences in their 
SFRs and Tl{U) values may be explained simpl y by differ- 
ent am ounts of gas available for star formation ( Kennicutt 



1998a 



drops below relation (6) somewhere in the magnitude range et al. 1995 



and NGC 3256 ( Zepf et al. 1999 ), but also in starbursts 
t hat are not directly r elated to merger events like e.g . M82 
( [O'Connell et al. 1995| ) and NGC 5253 ( [Gorjian 199"^ ). Such 
galaxies contain large numbers of highly luminou s young 
clusters, emitting up to 15% of the blue or U V light ( Meurer 
epf et al. 199^ and mark an extension of the 



Ssfr-Tl(C/) relation (LR2000). At the other extreme of the 
relation are galajcies like IC 1613 with exceedingly low (but 
non-vanish ing) star formation rates and very few star clus- 
ters at all ( [Wyder, Hodge fc Cole 2000| ). 

It thus appears that galaxies form YMCs whenever the 
SFR is high enough. This may have important implications 
for the understanding of how massive clusters (aka globu- 
lar clusters) formed in the early Universe. At earlier epochs, 
the general level of star formating activity was presumably 
higher because of the larger amounts of gas available, and 
globular clusters may have formed quite naturally. The high- 
est levels of star forming activity may have existed at the 
centres of rich galaxy clusters, giving rise to the rich popu- 
lations of globular clusters seen around many cD galaxies. 

The key to understanding massive cluster formation 
may lie in the properties of the interstellar medium in their 
parent galaxies. In the Milky Way, star clusters form in the 
cores of highly fra gmented giant molecular clouds (GMCs) 
(Lada et al. 1997), but at a very low efficiency averaged 
over the entire GMC. The largest GMCs in the Milky Way 
have masses of about 5 x lO^M© with a sharp cut-off above 



this limit (McKee 1999). It is tempting to speculate that 
this upper cut-off may be related to the upper limit of the 
open cluster luminosity function discussed in the previous 
section. If massive clusters form with similar low efficien- 
cies, their parent clouds must be much larger than Milky 
Way GMCs and it has been suggested that globular clus- 
ters formed from suyergia nt molecular clouds (SGM Cs) with 
masses of 10* - IO^Mq ( [Harris fc Pudritz 1994| ). In M51 
and NGC 5236, high-resolution CO studies have in fact re- 
vealed "Giant Mo lecular Associations" (CM As) with masses 



of 10^ - 10* Mr (Vogel et al. 1988; Rand fc Kulkarni 199C 



Rand et al. 1999), which might be iden tifiab le as the birth- 
sites of YMCs. Recently, Wilson et al. (2000) reported even 
more massive GMAs with masses of 3 — 6 x 10* M© in the 
YMC-rich "Antennae" galaxies. Such GMAs may assemble 
more easily in galaxies with high g as densities and, conse- 
quently, high star formation rates (Kennicutt 1998a). The 
connection between SFR and YMC richness could then be 
understood as resulting from an underlying dependence of 
both on the gas density. 



7 CONCLUSIONS 

Young Massive Star Clusters (YMCs) have been identified 
in the nearby face-on Sbc-type spiral M51. The richness of 
the YMC population in M51, as measured by the specific U- 
band luminosity Tl{U) is comparable to other cluster-rich 
spiral galaxies, compatible with the high star formation rate 
in M 51 and the Esfr - Tl{U) relation (Larsen fc Richtler 



LR2000). M51 itself is also very gas-rich (see e.g. ta- 
ble 1 in Kennicutt 1998b), consistent with its relatively high 
Esfr and Tl{U) values. 



200C). The high level of star formation activity may in this 
particular case be partly stimulated by interaction with the 
nearby companion NGC 5195, but interactions do not in 
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general seem to be a necessary requirement for YMC forma- 
tion. 

Within order-of magnitude estimates, the number of 
YMCs in M51 is compatible with extrapolation of a power- 
law luminosity function with exponent a ~ — 2 from a popu- 
lation of low-mass open clusters similar to that in the Milky 
Way. Thus, YMCs may simply represent a continuation of 
the normal open cluster luminosity function, extending to 
different upper limits in different galaxies. Therefore YMCs 
plausibly form in much the same way as lower-mass clusters 
and it appears likely that these same basic mechanisms also 
applied to the formation of globular clusters in the early 
Universe. 
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Table 3. Data for clusters in M51. My, U-B and B-V have 
been corrected for a foreground extinction of As = 0.150 mag. 
Clusters marked with an asterisk (*) were included on the HST 
F439W and F555W pointings. 
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Table pi (continued) 



ID 


RA(2000.0) 


DEC(2000.0) 


Mv 


U-B 


B-V 


984 


1 3-3n-0.5 03 


47:12:33.1 


— 11.12 


—0.88 


0.08 


QS7 
yo ( 


1 '^■'in-n^ OA 

10.OU.U0.U4: 


47-1 n-1 fi 


Q ^9 


— l.UO 


U.U I 


993 


13:30:05.19 


47:10:50.5 


-8.85 


-0.22 


0.27 


1013 


13:30:05.70 


47:09:49.9 


-9.90 


-0.32 


0.27 


1019 


13:30:05.92 


47:11:40.5 


-9.82 


-0.62 


0.35 


1031 


13:30:06.55 


47:11:01.1 


-9.51 


-0.75 


0.21 


1033 


13:30:06.67 


47:14:20.7 


-11.57 


-0.86 


0.08 


1049 


13:30:07.40 


47:11:26.5 


-10.12 


-1.01 


0.03 


1066 


13:30:08.47 


47:13:59.7 


-9.34 


-0.22 


0.36 
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